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ABSTRACT 

A micromechanical modeling procedure is developed to predict the viscoelastic properties of a graphite nanoplatelet/epoxy 
composite as a function of volume fraction and nanoplatelet diameter. The predicted storage and loss moduli from the model 
are compared to measured values from the same material using Dynamical Mechanical Analysis, nanoindentation, and tensile 
tests. In most cases, the model and experiments indicate that for increasing volume fractions of nanoplatelets, both the 
storage and loss moduli increase. Also, in most cases, the model and experiments indicate that as the nanoplatelet diameter 
is increased, the storage and loss moduli decrease and increase, respectively. 

INTRODUCTION 

In recent years, nanostructured materials have spurred considerable interest in the materials research community partly 
because of their potential for large gains in mechanical properties relative to current materials used for aerospace applications. 
This is particularly important in the design and development of Unmanned Aerial Vehicles (UAVs), in which the primary 
requirements are long-duration, high-altitude flights. Graphite nanoplatelet/polymer composites have been recently developed 
for use as a low-cost, lightweight materials with properties potentially superior to those of pure polymers. In order to facilitate 
the development of these materials for UAVs, modeling and characterization procedures must be developed over the nano- to 
macro-length scales. 

Recently, experiments were conducted to determine the viscoelastic properties of a series of graphite nanoplatelet/epoxy 
composite specimens that were fabricated at Michigan State University [1]. The specimens contained varying amounts of 
nanoplatelet volume fractions, sizes, and preparations. The experimental tests included nanoindentation, Dynamical 
Mechanical Analysis (DMA), and tensile testing. Since these tests were conducted over a wide range of length scales, the 
viscoelastic properties for each test must be interpreted based on the relative scale of the test method and scale of the 
heterogeneity of the material. However, the interpretation and comparison of the experimental data is difficult without the 
knowledge of the expected viscoelastic properties. 

In the current paper, a modeling procedure has been developed to facilitate the interpretation of the experimental data. A brief 
summary of the material and test methods is followed by a detailed explanation of the modeling procedure. The predicted 
viscoelastic properties from the modeling are compared to the experimentally obtained properties. 

MATERIALS 

In this study, the viscoelastic properties of four graphite nanoplatelet/epoxy composite materials were examined. The four 
materials had two different average sizes of graphite nanoplatelets, 1 and 15 pm diameters, for two different nanoplatelet 
volume fractions, 0.5% and 3.0%. The test specimens were fabricated at Michigan State University, and further details on the 
material fabrication and characterization can be found elsewhere [2]. Fig. 1 shows an Atomic Force Microscope (AFM) image 
of the 1 pm-reinforced composite with a nanoplatelet volume fraction of 3%. From Fig. 1 it is clear that the nanoplatelets are 
well dispersed and randomly oriented in the epoxy matrix. 




TEST METHODS 

Three test methods were used to determine the mechanical properties of the four composite materials, Dynamical Mechanical 
Analysis (DMA), quasi-static tensile testing, and nanoindentation. In DMA testing, an oscillatory force is applied onto a bulk 
sample and the resultant dynamic modulus is measured. These tests were standard tests based on well-established 
procedures [3, 4]. Quasi-static tensile tests were also performed on bulk samples. Using a standard tensile-testing machine 
and strain-measuring equipment, the Young’s modulus and Poisson’s ratio were determined for the four composite materials. 
Nanoindentation tests were performed on samples of the four materials using a previously established procedure [5]. In the 
nanoindentation procedure, a nanometer-sized tip was indented up to 2.0 pm into the specimen. A relatively small oscillatory 
displacement was superimposed onto the overall tip displacement, and the resulting dynamic modulus was measured. Further 
details on the application of these test methods to the four composite materials considered in the study may be found 
elsewhere [1]. 
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Fig. 2. Micromechanical modeling of graphite nanoplatelet/epoxy composite 



MICROMECHANICS MODEL 


The viscoelastic properties of the four composite materials were predicted by using a micromechanics model that incorporated 
the viscoelastic properties of the graphite and epoxy constituents. The micromechanics-based Mori-Tanaka method [6, 7] was 
used in conjunction with the elastic-viscoelastic correspondence principle [8]. To implement this model for the graphite 
nanoplatelet/epoxy composite, it was assumed that the nanoplatelets were randomly-oriented oblate-spheroids that were 
perfectly bonded to the surrounding matrix (Fig. 2). In this case, the composite complex modulus is 
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where c p and c m are the nanoplatelet and matrix volume fractions, respectively, C* and are the complex stiffness tensors 
of the nanoplatelets and matrix, respectively, I is the identity tensor, the angle-brackets indicate a nanoplatelet-orientation 
average, and A* is the complex dilute strain-concentration tensor of the nanoplatelets. The complex quantities in Eqn. (1) 

and throughout this paper can be expressed in terms of storage (time-independent) and loss (time-dependent) components. 
For example, the complex stiffness tensor of the composite is 

C* = C' + /C" (2) 

where C' and C" are the composite storage and loss stiffness tensors, respectively. The complex dilute strain-concentration 
tensor in Eqn. (1) is 
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where S* is the complex Eshelby tensor. The complex Eshelby tensor is determined by applying the elastic-viscoelastic 
correspondence principle to the Eshelby tensor for elastic composites [9, 1 0], which yields 
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where v* is the complex Poisson’s ratio; a x = a 2 > a 3 are the principle half-axes of the oblate spheroid inclusion (Fig. 2); 
S* jkl = S* jikl =S * jlk ; all other non-zero components are obtained by cyclic permutation of (1,2,3); the components which cannot be 
obtained by the cyclic permutation are zero; and 
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For three-dimensional randomly-oriented nanoplatelets, the orientation average of a tensor, T, is 
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where ij,k,l = 1,2,3; the indicial summation convention is used; 8^ is the Kronecker delta; and 
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Therefore, from Eqns. (6) and (7), (t) is isotropic. 


The properties of the constituent materials used in the modeling are shown in Table 1 . For this study, it was assumed that the 
graphite nanoplatelets were linear-elastic, with the classic transverse-isotropic graphite mechanical properties being used [11], 
where the x 3 axis is perpendicular to the graphene plane. The storage and loss moduli, E' and E " , respectively, of the epoxy 
material were determined via DMA testing [1]. It was assumed that the storage shear modulus, G ' , was equal to the static 
shear modulus, which was computed using the Young’s modulus and Poisson’s ratio determined via the quasi-static tensile 
testing [1]. The loss shear modulus, G " , was determined using [12] 
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With these assumptions, it was possible to uniquely determine all of the components of the complex stiffness tensors of the 
nanoplatelets and polymer using the elastic-viscoelastic correspondence principle. 


Table 1 . Viscoelastic properties of constituent materials in GPa 
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Epoxy 
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For the complex Eshelby tensor, the complex properties of the epoxy and the shape of the nanoplatelets were used in Eqns. 
(4) and (5). From AFM images of the material, such as that shown in Fig. 1 , it was determined that the graphite nanoplatelets 
had a thickness of about 0.5 pm for the 1 pm nanoplatelets and 1 pm for the 15 pm nanoplatelets [1]. 

It is noted that even though it has been previously shown that the mechanical properties of composites with nanometer-sized 
reinforcement must be modeled using an equivalent-continuum modeling approach [13, 14], the smallest dimension of the 
nanoplatelet is 0.5 pm, which is large enough to be modeled directly using micromechanics, as demonstrated previously [15]. 

RESULTS 

The storage moduli of the four composite materials measured in the DMA and nanoindentation experiments and predicted with 
the model are shown in Fig. 3. In addition, the Young’s moduli of three of the composite materials as measured via the tensile 
tests are shown in Fig. 3. The loss moduli of the composite materials measured through DMA and nanoindentation 




experiments and predicted by the model are shown in Fig. 4. The error bars shown for the experimental data in Figs. 3 and 4 
are the standard error as measured through repeat tests. 

It is evident in Fig. 3 that there is reasonable agreement between all of the experimental tests. Except for the 0.5% volume 
fraction with 1pm nanoplatelets, the DMA values of storage modulus are higher than the nanoindentation values for each 
material. There is no clear trend for the three tensile test values. Further analysis of the data leads to two important points. 
First, for all three experimental sets of data, the storage modulus (or Young’s modulus for the tensile tests) increases with an 
increasing nanoplatelet volume fraction, as expected. The model confirms this trend. Second, for the nanoindentation and 
tensile test data, the storage (or Young’s) modulus increases with the decreasing nanoplatelet size. The model confirms this 
trend, even though the DMA data does not follow this trend. Therefore, decreasing the diameter (while decreasing the 
thickness by a lesser amount), as done in the modeling, results in increasing storage (or Young’s) moduli. 

In Fig. 4, it is clear that the loss moduli measured with nanoindentation are at least 50% higher than those measured with the 
DMA for the four composite materials, and that the model shows much better agreement with the DMA data than for the 
nanoindentation data. The latter trend is expected because the model used the epoxy viscoelastic properties obtained from 
DMA tests. While the experimental data show an increase in loss modulus with increased nanoplatelet volume fraction, the 
model does not exhibit this trend with the 1 pm nanoplatelets. Furthermore, both experiments and model indicate an increase 
in loss modulus with increased nanoplatelet diameter. 
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Fig. 3. Storage moduli for graphite nanoplatelet/epoxy composites 
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Fig. 4. Loss moduli for graphite nanoplatelet/epoxy composites 


SUMMARY 

In this study, a micromechanical modeling procedure was developed to predict the viscoelastic properties of graphite 
nanoplatelet/epoxy composites as a function of volume fraction and nanoplatelet diameter. The objective of the model was to 
compare the predicted with the experimentally measured viscoelastic properties of this material from DMA, nanoindentation, 
and tensile tests, and thus facilitate the interpretation of the experimental data. 

The results of the modeling and all three sets of experiments indicate that for increasing volume fraction of nanoplatelets, the 
storage modulus increases. Even though an increase in nanoplatelet volume fraction leads to higher values of loss modulus in 
the experiments, the model does not show this trend. Except for the DMA data, the tests and model indicate that the storage 
modulus increases with decreasing nanoplatelet. For the same decrease in nanoplatelet diameter, the modeling and 
experiments both show a decrease in loss modulus of the composite. 
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